Hastatic Order in URu2Si2 
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Abstract 

The development of collective long-range order via phase transitions occurs by the 
spontaneous breaking of fundamental symmetries. Magnetism is a consequence of bro- 
ken time-reversal symmetry while superfluidity results from broken gauge invariance. 
The broken symmetry that develops below 17. 5K in the heavy fermion compound 
URu2Si2 has long eluded such identification. Here we show that the recent observation 
of Ising quasiparticles in URu2Si2 results from a spinor order parameter that breaks 
double time-reversal symmetry, mixing states of integer and half-integer spin. Such 
"hastatic order" hybridizes conduction electrons with Ising 5/ 2 states of the uranium 
atoms to produce Ising quasiparticles; it accounts for the large entropy of conden- 
sation and the magnetic anomaly observed in torque magnetometry. Hastatic order 
predicts a tiny transverse moment in the conduction sea, a collosal Ising anisotropy 
in the nonlinear susceptibility anomaly and a resonant energy-dependent nematicity 
in the tunneling density of states. 



The hidden order (HO) that develops below 17.5K in the heavy fermion compound 
URu2Si2 is particularly notable, having eluded identification for twenty five years [1-13]. Re- 
cent spectroscopic [14-18], magnetometric[19] and high-field measurements [20, 21] suggest 
that the HO is connected with the formation of an itinerant heavy electron fluid, a conse- 
quence of quasiparticle hybridization between localized, spin-orbit coupled f-moments and 
mobile conduction electrons. Though the development of hybridization at low temperatures 
is usually associated with a crossover, in URu 2 Si 2 both optical [22] and tunnelling [15-17] 
probes suggest that it develops abruptly at the HO transition, leading to proposals[10, 11] 
that the hybridization is an order parameter. 

High temperature bulk susceptibility measurements on URu2Si2 show that the local 5f 
moments embedded in the conduction Ising in nature [1, 26], while quantum oscil- 

lation (QO) experiments deep within the HO phase reveal that the quasiparticles exhibit a 
giant Ising anisotropy[21, 23, 24]. The Zeeman splitting depends solely on the c-axis com- 
ponent of the magnetic field, AE = g(9)fi B B[2A], with a g-factor g(9) = gcosO, where 6 
is the angle relative to the c-axis. The g-factor anisotropy exceeds 30, corresponding to 
an anisotropy of the Pauli susceptibility in excess of 900; this anisotropy is also observed 
in the angle-dependence of the Pauli-limited upper critical field of the superconducting 
state[23, 24], showing that the Ising quasiparticles pair to form a heavy fermion supercon- 
ductor. This giant anisotropy suggests that the f-moment is transferred to the mobile quasi- 
particles through hybridization [28]. In the tetragonal crystalline environment of URu 2 Si 2 
such Ising anisotropy is only natural in integer spin configurations^, 27], thus the most likely 
valence of the U ions with an integer spin is a 5/ 2 configuration. Moreover, the observation 
of paired Ising quasiparticles in a superconductor with T c ~ 1.5K indicates that this 5/ 2 
configuration is degenerate to within an energy resolution of gjj,BH c ~ 5K. In URu 2 Si 2 , 
tetragonal symmetry protects a magnetic non-Kramers T 5 doublet, the candidate origin of 
the Ising quasiparticles [4, 29]. 

The quasiparticle hybridization of half-integer spin conduction electrons with an integer 
spin doublet in URu 2 Si 2 has profound implications for hidden order; such mixing can not 
occur without broken double time-reversal symmetry. Time-reversal, O, is an anti-unitary 
quantum operator with no associated quantum number [30]. However double-reversal © 2 , 
equivalent to a 2n rotation, forms a unitary operator with an associated quantum number, 
the "Kramers index" if [30]. The Kramers index, K = (— 1) 2J of a quantum state of to- 
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tal angular momentum J defines the phase factor acquired by its wavefunction after two 
successive time-reversals, Q 2 \tp) = K\ip) = \vp 27T ). An integer spin state \ct) is unchanged 
by a 2tt rotation, so \a 2n ) = +\oi) and K — 1. However, conduction electrons with half- 
integer spin states, |ker), where k is momentum and a is the spin component, change sign, 
Ika 2 ") = -|k(j), soK = -1. 

While conventional magnetism breaks time-reversal symmetry, it is invariant under O 2 
so the Kramers index is conserved. However in URu2Si2, the hybridization between in- 
teger and half-integer spin states requires a quasiparticle mixing term of the form T-L = 
(\ka)V aa (k) (a\ + H.c) in the low energy fixed point Hamiltonian. After two successive 
time-reversals 

\ka)V aa (k)(a\ -+ |k^)Ok)(« 2 1 = -|k<7>y ff 2 :(k)(«|. (1) 

Since the microscopic Hamiltonian is time- reversal invariant, it follows that V aa (k) = 
— V 2 £(k); the hybridization thus breaks time- reversal symmetry in a fundamentally new 
way, forming an order parameter that, like a spinor, reverses under 2ir rotations. The re- 
sulting "hastatic (Latin: spear) order", is a state of matter that breaks both single and 
double time-reversal symmetry and is thus distinct from conventional magnetism. 

Indirect support for time-reversal symmetry-breaking in the hidden order phase of 
URu2Si2 is provided by recent magnetometry measurements that indicate the development 
of an anisotropic basal-plane spin susceptibility, Xxy, at the hidden order transition [19]. Xxy 
is a conduction electron response to scattering off the hidden order (c.f. Fig. 2.), lead- 
ing to a scattering matrix of the form i(k) = (<j x + a y )d(k), where d(k) is the scattering 
amplitude. t(k) has been linked to a spin nematic state [12], under the special condition 
that d(— k) = — d*(k) to avoid time-reversal symmetry breaking. However, if the scatter- 
ing process involves resonant hybridization in the f-channel, then d(k) is associated with 
resonant scattering off the f-state, a process with a real, even parity scattering amplitude, 
d(k) = d(— k). In this case, the observed t-matrix is necessarily odd under time-reversal in 
the hidden order phase. 

Hybridization in heavy fermion compounds is usually driven by valence fluctuations mix- 
ing a ground-state Kramers doublet and an excited singlet (cf. Fig. la). In this case, the 
hybridization amplitude is a scalar that develops via a crossover, leading to mobile heavy 
fermions. However valence fluctuations from a 5f 2 ground-state create excited states with 
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an odd number of electrons and hence a Kramers degeneracy (cf. Fig. lb). Then the quasi- 
particle hybridization has two components, ^/ a , that determine the mixing of the excited 
Kramers doublet into the ground-state. These two amplitudes form a spinor defining the 
"hastatic" order parameter 

-(!)• (2) 

The presence of distinct up/down hybridization components indicates that \& carries the 
global spin quantum number; its development must now break double time-reversal and 
spin rotational invariance via a phase transition. 

Under pressure, URu2Si2 undergoes a first-order phase transition from the hidden order 
(HO) state to an antiferromagnet (AFM) [31]. These two states are remarkably close in 
energy and share many key features[20, 32, 33] including common Fermi surface pockets; 
this motivated the recent proposal that despite the first order transition separating the two 
phases, they are linked by "adiabatic continuity," [32] corresponding to a notional rotation 
of the HO in internal parameter space [5, 34]. In the magnetic phase, this spinor points 
along the c-axis 

^ ~ Q , #b ~ ^ j (3) 

corresponding to time-reversed configurations on alternating layers A and B, leading to a 
large staggered Ising moment. For the HO state, the spinor points in the basal plane 

1 /e-^/A 1 /-e-^/A 

where again, = O^a and it is protected from developing a large moment by the pure 
Ising character of the 5/ 2 ground-state. 

Hastatic order permits a direct realization of the adiabatic continuity between the HO 
and AFM in terms of a single Landau functional for the free energy 

f[T, P, B z ] = [a{T c — T) — Vz Bl\ |*| 2 + /3|*| 4 - 7 (*V^) 2 (5) 

where 7 = S(P — P c ) is a pressure-tuned anisotropy term. The unique feature of the theory 
is that the non-Kramers doublet has Ising character, and only couples to the z-component 
of the magnetic field B z = B cos 9. The resulting Ising splitting of the non-Kramer's doublet 
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suppresses the Kondo effect, giving rise to the B 2 term in the quadratic coefficient, where 
the coefficient r\ z is of order 1/T^ Q . [35] The phase diagram predicted by this free energy 
is shown in Fig. 1 (c). When P < P C {T), the vector ^aty = \^\ 2 (n x , n y , 0) lies in the 
basal plane, resulting in hastatic order. At P = P c , there is a first order "spin-flop" into an 
magnetic state where = \^\ 2 (0, 0, ±1) lies along the c-axis. 

Adiabatic continuity provides a natural interpretation of the soft longitudinal spin fluctu- 
ations observed to develop in the HO state [36] as an incipient Goldstone excitation between 
the two phases [34]. In the HO state, rotations between hastatic and AFM order will lead 
to a gapped Ising collective mode which we identify with the longitudinal spin fluctuations 
observed in inelastic neutron scattering [36]. At the first order line, P = P c , the quartic 
anisotropy term vanishes; we predict that the gap to longitudinal spin fluctuations will soften 
according to A oc a/t]^P ~ \^\^/P c (T) — P (see Supplementary material). Experimental 
observation of this feature would provide confirmation of the common origin of the hidden 
and AFM order. 

Another prediction of the phenomenological theory is the development of a nonlinear sus- 
ceptibility anomaly with a collosal Ising anisotropy. From the Landau theory [35], the jump 
in the specific heat AC, the susceptibility and nonlinear susceptibility anomalies dxi/dT and 
A% 3 obey the relation ^A%3 = 12 (^-) 2 , where ^ = — ^r] z cos 2 9, so that % 3 oc cos 4 #. 
A large anomaly in the c-axis nonlinear susceptibility of URu 2 Si 2 has been observed at 
T c [25, 26], but its Ising anisotropy has never been quantified. The predicted development 
of a collosal Ising anisotropy in the zero-field non-linear susceptibility at the hidden order 
transition is another consequence of hastatic order. 

We now present a model that relates hastatic order to the valence fluctuations in URu2Si2, 
based on a two-channel Anderson lattice model. The uranium ground-state is a 5f 2 Ising 
T 5 doublet[4], |±) = a\ ± 3} + b\ =F 1), written in terms of J = 5/2 f-electrons in the three 
tetragonal orbitals Tf and T 6 



The lowest lying excited state is most likely the 5/ 3 (J = 9/2) state, but for simplicity 
here we take it to be the symmetry-equivalent 5/ 1 state. Valence fluctuations from the 
ground state (5f^ T 5 ) to the excited state (5f : r+ ) occur in two orthogonal conduction 



1+) = (<V^ + <t/r- t )l°> 




(6) 
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channels, [37, 38] T 7 and T 6 . This allows us to read off the hybridization matrix elements of 
the Anderson model 

H VF (j) = ^ 6 4 6± (i)|r+±)(r 5 ± | + v 7 cl T ( 3 )\rlT)(r 5 ± | + h.c. (7) 

where ± denotes the "up" and "down" states of the coupled Kramers and non-Kramers 
doublets. The field c\. a {j) = J2k ^rOO c^ T e~ lk ' Kj creates a conduction electron at site 

L J tJT 

j with spin a in a Wannier orbital with symmetry Y G {6,7}, while V e and V 7 are the 
corresponding hybridization strengths. The full model is then written 

H = E e * c l c *° + E 1 H vfU) + H a (j)] (8) 

kcr j 

while H a (j) = AEJ2± |r 7 d=, j)(T 7 ±,j\ is the atomic Hamiltonian. 

Hastatic order is revealed by factorizing the Hubbard operators, |r7(j}(r 5 a| = ^>\xa- 
Here \T 5 a) = xLlty is the non-Kramers doublet, represented by the pseudo-fermion x a j 
while is a slave boson[39] representing the excited f 1 doublet \T^a) = ^{Q). Hastatic 
order is the condensation of this boson ^Ixa — > (&o)Xon generating a hybridization be- 
tween the conduction electrons and the Ising 5/ 2 state while also breaking double time 
reversal. The T 5 doublet has both magnetic and quadrupolar moments represented by 
X^&X = {Ox 2 - y 2 ,Oxy, mZ ), where m z is the Ising magnetic moment and O x 2_ y 2 and O xy 
are quadrupole moments. The tensor product Q a/ 3 = ^o^^ describes the development of 
composite order between the non-Kramers doublet and the spin density of conduction elec- 
trons. Composite order has been considered by several earlier authors in the context of two 
channel Kondo lattices[37, 42, 43] in which the valence fluctuations have been integrated 
out. However, by factorizing the composite order in terms of the spinor \I/ a , we are able to 
directly understand the development of coherent Ising quasiparticles and the broken double 
time-reversal. 

Using this factorization, we can rewrite the valence fluctuation term as, 

HvfU) = E c L^(k,i)x„0>- ik ^ + H.c. (9) 



. * t . 

where V(k, j) = V 6 $ r J(k)B\ + V^r- ( k ) B\o x , B 3 = \ .In the ordered state, 

,: tf, ' 
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Bj = (Bj) is replaced by its expectation value, so that in the HO state 

/ e i(Q-Rj+*)/2 o \ 

< B 1> = I*I YOR+il/2 hl% (io) 

I e -«(Q-Rj+<^)/2 / 

with magnitude \^\. The internal angle, 0, rotates -Bj within the basal plane. 

As the HO and AFM appear to share a single commensurate wavevector, Q = 
(0, 0, ^f)[20, 32, 33], we use this wavevector here. It is convenient to absorb the unitary 
matrix, Uj into the pseudo-fermion, Xj — UjXj- This gauge transformation transfers the 
charge from the slave boson to the pseudo-fermion, making it a charged quasiparticle. In 
this gauge, one channel (T 6 ) is uniform, while the other (Tj) is staggered 

Hvf = c k V e(k)xk + ctV 7 (k)xk +Q + h.c. (11) 

k 

where the hybridization form factors VV(k) = V7$7^(k)<7i and Ve(k) = Vq^^O^)- 

There are two general aspects of this condensation that deserve special comment. First, 
the two-channel Anderson impurity model is known to possess a non-Fermi liquid ground- 
state with an entanglement entropy of |/c B ln2[40]. The development of hastatic order 
liberates this zero-point entropy, accounting naturally for the large entropy of condensation. 
As a slave boson, \1> carries both the charge e of the electrons and the local gauge charge 
Qj = ^j^j + XjXj °f constrained valence fluctuations, its condensation gives a mass to their 
relative phase via the Higg's mechanism[41] . But as a Schwinger boson, the condensation 
of \1/ this process breaks the SU (2) spin symmetry. In this way the hastatic boson can be 
regarded as a magnetic analog of the Higg's boson. 

One of the key elements of the hastatic theory is the formation of mobile Ising quasipar- 
ticles, and the observed Ising anisotropy enables us to set some of the parameters of the 
theory. The full anisotropic g-factor is a combination of f-electron and conduction electron 
components given by g(9) ~ g* cos# + g c (^=r) where g c = 2 and the factor Tk/D derives 
from the small conduction electron admixture in the quasiparticles. Experimentally [21], the 
g-factor anisotropy is in excess of 30, i.e. g z /g± ~ D/T K > 30, which enables us to phe- 
nomenologically set a lower bound on D/T K in our model. The g-factor is defined in terms 
of the Zeeman splitting of the heavy fermion dispersion, AE kr) = \E kri ^ — E kr) ±\ = g^i^^BB, 
where rj G [1, 4] is a band-index. Fig. 3(a) shows the Fermi surface averaged g-factor, defined 
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by 



9(0) 



Eta, WW- 5 **?) 



(12) 



calculated within the mean-field hastatic model, as a function of field angle to the c-axis, 
choosing the lower-bound estimate D/T K ~ 30. 

Another key aspect of the hastatic picture is that there must be time-reversal symmetry 
breaking in both the HO and AFM phases, manifested by a staggered moment; in the AFM 
this leads to a large c-axis f-electron moment, but in the HO it becomes a small transverse 
moment carried by conduction electrons, m c = —gfiBTrcrG c (k, k + Q), where Q c is the con- 
duction electron Green's function[35]. The small magnitude of the induced moment is a 
consequence of the Clogston-Anderson compensation theorem that states that changes in 
the conduction electron magnetization are set by the same ratio Tk/D that determines the 
g-factor anisotropy[44]. There will also be a small mixed- valent contribution from the excited 
Kramers doublet, rh\ oc (Waty). The angle of the moments in the plane is controlled by 
the internal hastatic angle, 0. Fig. (3b) shows the temperature dependence of the magnetic 
moment calculated for a case where D/Tk ~ 30, for which m_i_(0) =0.015^, which is an 
upper bound for the predicted conduction electron moment. Neutron scattering measure- 
ments on URu2Si2 have placed bounds on the c-axis magnetization of the f-electrons using 
a momentum transfer Q in the basal plane. Detection of an m_|_(0) carried by conduction 
electrons, with a small scattering form-factor, will require high-resolution measurements 
with a c-axis momentum transfer. We note that there have been reports from /iSR and 
NMR measurements [45, 46] of very small intrinsic basal plane fields in URu2Si2, which are 
consistent with this theory. 

Although the conduction electrons develop a magnetic moment, the non-Kramers 51 2 
state carries no dipolar or quadrupolar moment, (x^&X) = 0- I 11 the microscopic model the 
quadrupolar moments vanish because of the d-wave form factor between the T 6 and T 7 - 
channels[35]. The absence of a charge quadrupole implies there will be no lattice distortion 
associated with hastatic order. By contrast, hastatic order does induce a weak broken 
tetragonal symmetry in the spin channel. In the HO state, the inter-channel components 
of the hastatic t-matrix, V§V\ oc a x + <J y , break tetragonal symmetry in the spin channel, 
resulting in a nonzero spin susceptibility within the conduction fluid, 



X*y = -(gfi B ) 2 TTa x g c (k, k + QK£ c (k + Q, k) oc (TrV^) 2 



(13) 
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of a magnitude of order (^) , that onsets at the hidden order temperature as |\1/| 4 ~ 
(T c - Tf as shown in Fig. (3c). 

Hastatic order also manifests itself as an anisotropic hybridization gap, which vanishes 
along lines in momentum space, giving rise to a "V" -shaped density of states due to the 
partial gapping of the Fermi surface, as shown in Fig (4), which will be smeared out by 
disorder in the real material. The anisotropy of the hybridization also breaks tetragonal 
symmetry, giving rise to a energy dependent nematicity r](E) that peaks over a narrow 
energy window around the f-resonance. Some of this nematicity is present at the Fermi 
surface, accounting for the splitting seen in quantum oscillation frequencies [47] and cyclotron 
resonance experiments [48]. An ideal way to verify this prediction is via scanning tunneling 
spectroscopy (STS), where the measured differential conductivity, dl/dV oc A(eV, x) is 
proportional to the local density of states A(eV,x) at position x on the surface. A measure 
of the broken tetragonal symmetry is provided by the "nematicity" 

n(v) = §{x ' v) sign{xv \ . (U) 
((f) 2 -(f) 2 ) ' 

Here (x, y) is the co-ordinate relative to the center of the unit cell and the overbar denotes an 
average over the unit cell. The resonant scattering of the hastatic order causes this quantity 
to vary rapidly as a function of voltage, over an energy range of order the Kondo temperature 
T K . Fig. (4) shows the variation of the nematicity, calculated within our model of hastatic 
order. As one passes through the Kondo resonance, this quantity is found to change sign, 
and to peak at a maximum value of approximately 50%. 

So far we have focused on the mean field consequences of hastatic order. The recently 
observed softening of the commensurate longitudinal spin fluctuations at T c [49] suggests 
that hastatic order melts via phase fluctuations. While the hybridization spinor itself, (^), 
can only become nonzero below the phase transition at T c , we expect that its amplitude, 
(^ f ' 1 X f ), will persist to higher temperatures. As (Waty remains zero above T c , only the non- 
symmetry breaking (uniform, intrachannel) components of the hybridization can develop: 
V 6 Vg and V7V7 will emerge via a crossover at a higher temperature T* to create an incoherent 
Fermi liquid, consistent with the heavy mass inferred from thermodynamic and optical 
measurements [1, 18] and the development of Fano signatures in both scanning and point- 
contact tunneling spectroscopy [15-1 7]. The symmetry-breaking, interchannel components, 
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VeaiVj will always develop precisely at the HO transition. Another aspect of experiments 
that is not covered by our mean field description, is the observation of gapped, low-energy 
incommensurate fluctuations around a Q- vector Q* = (1 ± 0.4, 0, 0) in the hidden order 
phase, which appears to be a sign of an unfulfilled predisposition towards an incommensurate 
phase, probably driven by partial Fermi surface nesting. These effects lie beyond a mean- 
field description, but would emerge from the Gaussian fluctuations about the mean-field 
theory . 

Though we have discussed hastatic order in the context of URu 2 Si 2 , it should be a more 
a widespread phenomenon associated with hybridization in any f-electron material whose 
unfilled f-shell contains a geometrically-stabilized non-Kramers doublet As such, we expect 
realizations of hastatic order in other 5f uranium and 4f praseodymium intermetallics. 
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Figure Captions 

Fig. 1. (a) A normal Kondo effect occurs in ions with an odd number of f-electrons, where 
the ground state is guaranteed to be doubly degenerate by time-reversal symmetry 
(known as a Kramers doublet). Virtual valence fluctations to an excited singlet state 
are associated with a scalar hybridization, (b) In URu2Si2, quasiparticles inherit 
an Ising symmetry from a 5/ 2 non-Kramers doublet. Loss or gain of an electron 
necessarily leads to an excited Kramers doublet, and the development of a coherent 
hybridization is associated with a two-component spinor hybridization that carries 
a magnetic quantum number and must therefore develop at a phase transition, (c) 
Phase diagram for hastatic order, showing how tuning the parameter A oc (P — P c ). 
leads to a spin flop between hastatic order and Ising magnetic order. Inset: at the 1st 
order line, the longitudinal spin gap is predicted to vanish as A oc \JP C — P. (d) Polar 
plot showing the predicted cos 4 6 form of the non-linear susceptibility xz induced by 
hastatic order, where 6 is the angle of the field from the c-axis. 

Fig. 2. Phenomenological interpretation of the anomalous spin susceptibility in URu2Si2- (a) 
The anomalous spin susceptibility is given by conduction electrons scattering off the 
hidden order, sandwiched between a x and a y vertices, (b) The anomalous scattering 
t-matrix can be rewritten as a resonant scattering off the order parameter. 

Fig. 3. Magnetic response of hastatic order, (a) Polar plot of calculated g-factor, g{9) av- 
eraged over the Fermi surface, as a function of magnetic field angle 9 (see SOM for 
details), compared with results of Altarawneh et al. [21], overlaid in green, (b) As a 
consequence of the broken time-reversal symmetry, we predict a staggered conduction 
electron moment that onsets at the HO transition with a linear T c — T temperature 
dependence (staggering pattern shown in inset). The magnitude of this moment is gov- 
erned by Tk/D ~ .01/j.b/U. (c) We have calculated the tetragonal symmetry breaking 
component of the uniform susceptibility, Xxy(T). To compare our results to Okazaki 
et al [19] (overlaid as green squares), we have plotted the two-fold oscillation amplitude 
of the magnetic torque, A (in black), where A cos 20 = ^ = — cos 2<pXx y {T). A 
goes as (T c — T) 2 just below the HO transition. For details of our calculation, including 
parameter choices, please see the supporting online material [35]. 
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Fig. 4. Density of states and resonant nematicity predicted by theory. Upper panel: density of 
states as a function of energy predicted by model calculation (blue line), showing f and 
conduction electron components. Red line, voltage dependence of nematicity rj(V) in 
model calculation of scanning tunneling spectrum. Lower panels: spatial dependence 
of density of states for selected bias voltages in model calculation of scanning tunneling 
spectrum, showing the resonant character of the nematicity. 
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I. ABSENCE OF ISING SYMMETRY IN A KRAMERS DOUBLET WITH 
TETRAGONAL SYMMETRY. 



A central argument in the main paper is that the Ising symmetry of the quasiparticles in 
URu 2 Si 2 must derive from a non-Kramers doublet. This section establishes that a Kramers 
doublet in a tetragonal crystal lacks the selection rules required for an Ising symmetry 
without fine tuning. A tetragonal crystal field Hamiltonian contains terms of the form J± 
allowed by the four-fold symmetry in the basal plane. These terms cause a pure doublet 
| ± M) (where M is some J z G { — J, ■ ■ ■ J}), to mix with states | ± M') differing by four units 
of angular momentum, where M' = M — An, and n is any integer. A magnetic doublet in 
an tetragonal environment thus has the form 

|r+> =J2 a n\M-4n) 

n 

|r->=5> n |-M + 4n>. (1) 

n 

The coefficients a n may always be chosen to be real. Ising symmetry requires that the matrix 
elements 

(r + |J±|r-) = ^2a n a n ,(M -4n\J±\ — M + An') (2) 

n,n' 

vanish. In the absence of fine-tuning (a n = for all n ^ 0), this implies a selection 
rule (M - An\ J±\ - M + An') = 0. Such terms vanish if M - An ^ -M + An' ± 1, or 
M ^ 2{n — n') ± |. Since n and n' are integers, M must be an integer. For any half- 
integer M, corresponding to a Kramers doublet, the selection rule is absent and the ion will 
develop a generic basal plane moment. The fine-tuned case will produce an Ising Kramers 
doublet, but corresponds to the complete absence of tetragonal mixing, highly unlikely in a 
tetragonal environment. However, for integer M, corresponding to a non-Kramers doublet, 
the selection rule (M — An ± | J±\T — M + An') = causes every term in the above sum (2) 
to vanish so the transverse moment is necessarily zero, yielding perfect Ising symmetry. 

By contrast, in a hexagonal system like CeAl3, the crystal field Hamiltonian contains 
terms of the form J±. Such terms again mix a pure doublet | ± M) with terms | ± M'), 
where M' = M — 6n, integer n. However, for J < 7/2, there are no choices of M and 
M' that differ by 6, and thus there will be two Ising doublets for the Ce, J = 5/2 case: 
T 8 = | ± 5/2 and T 9 = | ± 3/2). Either of these Kramers doublets could undergo a single 
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channel Ising Kondo effect [9, 17], which will differ substantially from the two-channel Kondo 
physics associated with a non-Kramers doublet. 



II. ESTIMATE OF T K FOR URu 2 Si 2 

In our mean field theory, all Kondo behavior develops at the hidden order transition. 
Incorporating Gaussian fluctuations should suppress the hidden order phase transition, T H q, 
while allowing many of the signatures of heavy fermion physics, including the quenching of 
the spin entropy and the heavy mass to develop at a higher crossover scale, T K . While the 
coherence temperature estimated from the resitivity, T* 70K is much larger than the 
hidden order temperature, T HO = 17.5K, an entropic estimate of the Kondo temperature, 
S(Tk) = ^i21og2 gives an effective Kondo temperature not much larger than Tho- There 
is considerable uncertainty in the entropy associated with the development of hidden order, 
S(Tho), due to difficulties subtracting the phonon and other non-electronic contributions, 
leading to estimates ranging from .15i21og2[8] to .3/2 log 2 [7]. If we take a conservative 
estimate of S{T HO ) = .2i21og2, and the normal state 7 = 180mJ/molK 2 [7], S(T K ) = 
.2R log 2+ 7^T = log 2 yields T K = 27K, much lower than the coherence temperature 
seen in the resistivity. 

III. LANDAU THEORY FOR HASTATIC ORDER. 
A. Landau theory in zero field 

The most general Landau functional for the free energy density of a hastatic state with 
a spinorial order parameter f as a function of pressure and temperature is 



/[*] = a(T c - T)|Mf + m A ~ 7(*V 2 tf) 2 



(4) 



where 7 = S(P — P c ) is a pressure-tuned anisotropy term and 





where 9 is the disclination of Wa^ from the c-axis. Using this expression for 



/ = - a (T - T c )r 2 + /3r 4 - 7 r 4 cos 2 9. 



(6) 
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If P < P c , then 7 < and the minimum of the free energy occurs for 9 = tt/2, corresponding 
to the hidden order state ordered state. By contrast, if P > P c , then 7 > and the minimum 
of the free energy occurs at 9 = 0, n, corresponding to the antiferromagnet. The "spin flop" 
in 9 at P = P c corresponds to a first order phase transition between the hidden order and 
antiferromagnet (See Fig. 1) 

To study the soft modes of the hastatic order, we need to generalize the Landau theory to 
a time-dependent Landau Ginzburg theory for the action, with action S = J Ldtd 3 x, where 
the Lagrangian 



and p is the stiffness. Expanding \l/ around its equilibrium value \I/cb taking = for 
convenience and writing 



This gives rise to a change in ^a z ^ = x\^o\ 2 + 59(x)z\^q\ 2 corresponding to a fluctuation 
in the longitudinal magnetization. This rotation in ^ does not affect the first two isotropic 











P-Pc 



(3) 



FIG. 1. Global phase diagram predicted by Landau theory. 
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terms in f[^}. The variation in the action is then given by 

5S = p|* | 2 J2 l%)! 2 (V - ^ + ~ P )l^| 2 ) (8) 



The dispersion is therefore 



where 



a; 2 = {cqf + A 2 (9) 



, 25(P C -P) 



*o| 2 , (10) 

so that even though the phase transition at P = P c is first order, the gap for longitudinal 
spin fluctuations is 

Aoc \V Q \y/P c -P. 

Since dP c /dT c is finite, close to the transition, \JP C — P « ^/dP c /dT c (T — T c ), and A oc 
a/T — T c . Inelastic neutron scattering experiments can measure this gap a function of tem- 
perature at a fixed pressure where there is a finite temperature first order transition. The 
iron-doped compound, URu2_ :E Fe. c Si2 can provide an attractive alternative to hydrostatic 
pressure, as iron doping acts as uniform chemical pressure and tunes the hidden order state 
into the antiferromagnet[18]. 



B. Landau Theory in Magnetic Field: Nonlinear Susceptibility 

The origin of the large c-axis nonlinear susceptibility anomaly in URu2Si2 [10] has been a 
long-standing mystery. It has been understood phenomenologically within a Landau theory 
as a consequence of a large ^ 2 B 2 coupling of unknown origin[10, 11]. Here we show that 
the Landau theory of hastatic order will contain just such a term. 

While the conduction electrons couple isotropically to an applied field, the non-Kramers 
doublet linearly couples only to the z-component of the magnetic field, B z = B cos 9, which 
splits the doublet as it begins to suppress the Kondo effect. When we include the effect of 
the magnetic field in the Landau theory, we obtain: 

fm = [a(T c -T)-r ?ZJ B 2 -r ?±J Bi]^ 2 + /3|^| 4 + 7 (^V 2 ^)2 ! (n) 

where we shall show below that the coefficient of the ^> 2 B 2 Z term, r/ z goes as p/T^ Q , where 
p is the conduction electron density of states, while the coefficient of the ^> 2 B\ term, r]± is 
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of order p/D 2 , where D is the conduction electron bandwidth. Minimizing this functional 
with respect to we obtain 

/ = ~ [a(T c — T) — Vz (Bcos9) 2 - V± (B sin Of] 2 . (12) 

Following the arguments of [11], we can calculate the jump in the specific heat AC V and the 
susceptibility and non- linear susceptibility anomalies dxi/dT and A% 3 respectively, to find 

ACV a 2 



T HO 2/3 



(13) 



A X3 = | cos 2 + V± sin 2 fl) 2 « ^ cos 4 9 (15) 

where Xi an d X3 are the anomalous components of the susceptibility that develop at Tho- 
These results show that xi an d Xs wu l exhibit a giant Ising anisotropy. The thermodynamic 
relation 

—X3 = 12 — (16) 



dT 



still holds. 



1. Estimating rj z and rj± 

To illustrate this simple Landau theory, we will calculate fj in a simplified model: we will 
neglect the momentum dependence of both the f-level and the hybridization and take the 
hastatic order to be uniform. None of these assumptions qualitatively changes the results. 
The |\1/| 2 coefficient is calculated from the microscopic theory (see the next section) by 
expanding the action, S = — Trlog [l — ^(V^QoiV^)] in \I>, where To = (icu n — A — 
gffifBza-s)' 1 and Qo = (iu n — — g/2B ■ c?) _1 are the bare x an d conduction electron 
Green's functions (remember, x are the fermions representing the non- Kramers doublet). V 
represents the hybridization matrix elements, which are momentum-independent here, and 
proportional to the unit matrix. Note that while the conduction electrons are isotropic, the 
X's a re perfectly Ising. The coefficient of |\I/| 2 is then, 




^EE— — — ^ ( 17 ) 



ibj n — eko- iuj n — Ao- 



where eko- = e^ — g/2a • B and X a = A — gfp,faB z are the dispersions in field. Performing 
the Matsubara sum, we obtain 



v 2 y I deV{e) t&nh ^- tanh ^ = pV * y 



^ 18 

d 2(A a -e + g/23-B) 



where we approximated the conduction electron density of states as a constant, p within 
the bandwidth, 2D. Let us first calculate r]j_, quantizing the field along the x-direction and 
taking g = 2, 



-pV 2 I de 9 tanh 2T 



tanh 



IT 



D 



dB 2 



2(A - e + cjB a 



. (19) 



B, =0 



As the integrand is a function of e — aB, the integral is straightforward. And as D ^> A, T, 
the dominant term will be: 

D 



ij± =- p v 2 ' sech22T 



4T( 



A) 



+ 



tanh ^ — tanh ^ 



A) 2 



pV 2 



(20) 



-1? 



?7 2 will have three contributing terms: one purely from the conduction electrons that is 
one arising from cross terms between the conduction and f-electrons, and finally one solely 
from the f-electrons that dominates the other two. We shall focus on this last term, 




v - f D d 2 tanh ^ - tanh A " g y B ' 
-pV 1 y j I de 



L D dBl 2(\-e-g f p f aB z ) 
This integral cannot be done analytically at finite temperature, so we take T — Y 



(21) 



B z =0 



Vz 



D 



2(A C 



pV 2 



B z =0 



pV 2 

A 2 (D + A) 2 



pV 2 

1 HO 



(22) 



as A = Tho at zero temperature. So r)±/rj z = Using a conservative of T HQ /D ~ 1/30, 
we predict an anisotropy of about 900 in dxi/dT and nearly 10 6 in Ax3- However, in a 
realistic model, there will be f-electron contributions to n± involving fluctuations to excited 
crystal field states that may reduce the anisotropy somewhat. The important point here is 
that the anisotropies will be orders of magnitude larger than the single ion anisotropy in xi 
(approximately 3), and furthermore, that they will develop exclusively at the hidden order 
transition. 
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IV. TWO CHANNEL ANDERSON LATTICE MODEL FOR URU 2 SI 2 



Our model describes a lattice of U ions immersed in a conduction sea of electrons. We take 
the low energy configuration of each U ion to be a 5/ 2 T 5 non-Kramers doublet. All energies 
are measured relative to the energy of the isolated doublet. For simplicity, in our model we 
take the dominant valence fluctuation channel to be 5/ 2 # 5/ 1 . Particle-hole symmetry can 
be used to formulate the equivalent model with fluctuations into a 5/ 3 Kramers doublet. 
The full model is then written 

H = £ e * c l c ^ + £ (H VF (j) + H a (j)) (23) 

k<r j 

where cj^ creates a conduction electron of momentum k spin a, with energy ek, 

H VF (j) = v 6 4 6± (j)\r+±){r 5 ± | + v 7 4 7T (j)\r+T)(r 5 ± | + h.c. (24) 

describes the valence fluctuations between the T 5 doublet and an excited T 7 Kramers doublet 
while 

H a (j)=AE^\T 7 ±)(T 7 ±\ (25) 
± 

is the atomic Hamiltonian. 

In our model, we choose the T 7 to be the lowest excited state; in this case the the valence 
fluctuations are determined by the decomposition, 

|r 7 ±) = c#| T |r 5 ±> + /3^7 T |r 5 =F), (26) 

where tjr Ta = J2k c k/3 [^^\a(3 e~ ik ' Rj creates a conduction electron in a Wannier state with 
j = 5/2 and symmetry T (= r 6 , T 7 ) localized around the uranium atom at site j. 
For a plane wave, the form factors are given by [$rk] Qj g = ^^(k), where 

^(k) = Y 3m _« (k)(3m - f , \^\5/2m)&) 

In URu 2 Si 2 , the uranium atoms are located on a body centered tetragonal lattice (bet) at 
relative locations, Unn = (±a/2, ±a/2, ±c/2), while the f-electrons hybridize with conduc- 
tion electrons located around the silicon atoms, a^A? = (±a/2, ±a/2, ±z) where z = .371c 
is the height of the silicon atom above the U atom[3]. The form-factor is then, 

l^U= £ e-*-^yi^ NN ) (27) 

{Rjvjv,ajvjv} 
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Notice that this function has the following properties: [$rk+G] Q/3 = [$rk] Q/ 3 and [$rk+Q] a/3 = 

To cast the Hamiltonian as a field theory, we introduce a slave boson/slave fermion 
representation, |r 7 o")(r 5 a:| = b\xai where bl\£l) represents the occupation of 5/ 1 , and b a 
carries a positive charge and xjj^) represents the 5/ 2 state. The valence fluctuation term 
in the Hamiltonian then takes the form 

HvfiU) =Y,<Lv!&OL,j)XaU)e- <k - 1i > +H.c. (28) 
/$} . 

Introducing B* = (where we use W instead of the traditional b\ to reinforce 

\° *l 

that these are indeed order parameters), we can write the hybridization compactly, 

V (1) (k, j) = V^ u B) + Vr^-Bja,. (29) 

In the mean field theory, we replace V(k, j) — > (V(k,j)), replacing = \l>*. We 

consider the configuration 

e ^.w*) =mu > (30) 

where Uj is a unitary matrix, and Q = (0, 0, — ) is the common wavevector for hastatic order, 
chosen to match the antiferromagnet, and <fi = n/4 is chosen to match the susceptibility 
anisotropy. It is convenient to redefine Xj = UjXj- With this device, the spatial dependence 
of (Bj) is absorbed into the redefined f-electrons, so that BjXj = \^\Xj an d 

B)a lXi = \*m<riU})xj = |^|(n-a)e^- R ^„ n = -^(x + y). (31) 

In this gauge, the T 6 hybridization is uniform while the r 7 hybridization is staggered. We 
write V 6 (k) = |^|\/ 6 $r 6 , and V 7 (k) = |*|y 7 $r 7 (n • a). 

In the slave formulation, the atomic Hamiltonian is H a (j) = AEY2 a ^]a^ja- The intro- 
duction of slave bosons and fermions to represent the Hubbard operators requires a constraint 
to maintain one particle per site, (j2 a *5<r*j<r + E Q x] a Xja ~ lj ■ 

We take a simplified model of the conduction electron hopping, treating them as s-wave 
electrons located at the U site, hopping on a bet lattice with dispersion 

kr^d kyCL / \ 

ek = — St cos — - cos cos — \i. (32) 

Zi Zi Zi 



We do, however, want to capture the essential characteristics of the URu 2 Si 2 bandstructure 
- namely nesting between an electron Fermi surface about the zone center and a hole Fermi 
surface at Q[3]. In order to favor a staggered hybridization, and to match up with ARPES 
experiments suggesting a heavy f-band[2], we take the hole Fermi surface to be generated 
from a weakly dispersion \ band. This f-electron hopping will be naturally generated by 
hybridization fluctuations above T HO , effectively where (B^B) ^ while (B) = 0. A large N 
expansion of this problem would capture these fluctuation effects, but is overly complicated 
for this problem so we put this dispersion in by hand, e/k = — Stf cos %p cos ^ cos 
So to summarize, our mean-field Hamiltonian is, 

k k j j 

(33) 

We rewrite this Hamiltonian in matrix form 



H = J2 (ii+Q^LxL+q) 



+ ^[(A£ + A)^% T -A 





H 


,/3(k) 






( ek 





V 6 (k) 


V 7 (k) \ 


( ^ ^ 





e k+Q 


-V 7 (k) 


-V 6 (k) 


C k+Q 


V 6 f (k) 


-V{(k) 


A k 





Xk 


U(k) 


-v 6 f (k) 





^k+Q / 


\ Xk+Q ) 



(34) 



where we have suppressed spin indices, made the assumption that Aj = A is uniform, equiv- 
alent to enforcing the constraint on average, introduced Ak = A — e/k, and used the sim- 
plification that Q is half a reciprocal lattice vector, making V(k + Q) = — V(k), as shown 
above. 

In the absence of particle-hole symmetry, this Hamiltonian cannot be diagonalized ana- 
lytically, and must be done numerically, giving a set of four doubly degenerate bands, E^ v . 
The mean field free energy is, 

r 1 



F[b, A] 



log [1 + e~^ k "] +A/" S [2(AE + A)|^| 2 - A] , 



(35) 



k,r? 



where (5 = (k B T)~ l . The mean field parameters, and A are obtained by numerically 
finding a stationary point of the free energy that minimizes F with respect to \^/\ and 
maximizes it with respect to A. A plot of \^\ and A as a function of temperature is shown 
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in Figure 1 C, for the parameters used to calculate \xy{T) and m(T) for Figure in the main 
paper. The parameters are as follows: t = 12.5 meV is taken to match the magnitude 
of Xxy from the torque magnetometry data[6]. ///£ — —.075 gives the slight particle-hole 
asymmetry essential to reproduce the flattening of \ xy at low temperatures, and has also 
been adjusted so that /i + A = 0atT = for consistency with the dl/dV calculations; 
tf/t = —.025 gives a weak f-electron dispersion; the crystal field angle £ = .05 is taken 
to be small, as it is in CeRu2Si2[4]; Vq/Vj = 1 is arbitrary; and finally V 2 /AE = 2t is 
chosen to give 2|\E'| 2 15% mixed valency. The shape of the susceptibility curve is quite 
sensitive to the particle-hole asymmetry and degree of mixed valency, but not to the other 
parameters. The internal angle, <fi controls the susceptibility anisotropy and the orientation 
of the magnetic moments. = 7r/4 reproduces Xxy 7^ 0, Xxx = Xyy 

It is important to note that the hastatic state is promoted by a partial nesting of the 
normal state Fermi surface along Q = (100), leading to the partial gapping of the Fermi 



(A) 



;0,0,0^^<2jri,0) ^ 


(2^,iyr,0) 




/ (0,0,0) ' 


0,0,2* 




% j 


1 ? / 




r 




10 15 20 25 30 35 

(C) 



0.2 0.4 0.6 0.8 1.0 

T/T HO 



1.2 

/ 0.8 
/ 0.6 
/ 0.4 

0.2 


A(E) 






3 -2 -1 


1 3 E(meV) 4 



(36) 



FIG. 2. (A) Band structure of the hastatic order is shown in solid blue, while the bare conduction 
(red) and f (green) bands are dashed. (B) Mean field parameters |^| and A as a function of 
temperature. (C) Density of states in hastatic order, for the region close to the Fermi energy 
containing the hybridization gap. 
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surface seen in the density of states (Fig. 2C). However, even for perfect nesting, there will 
still be gapless excitations, as the hybridization gap, VqV 7 contains a node along the c-axis 
(due to a c-axis node in V7). In this sense (and this sense only), the hastatic state is similar 
to the Ikeda-Miyake model [16] of the failed Kondo insulators CeNiSn and CeRhSb, which 
also predicts a gapless state. 

Above Tjjo i n URu 2 Si 2 , there are magnetic fluctuations seen at two momentum trans- 
fers: at the commensurate Q = (1, 0, 0), which goes soft at Tho and at the incommensurate 
Q* = (1 ± .4, 0, 0), which remains at finite energy above both the HO and AFM phases[12]. 
While we have chosen a simplified bandstructure that captures the partial nesting at the 
commensurate Q, it is likely that a more realistic bandstructure will contain other ener- 
getically competitive nesting vectors that never develop due to the stronger nesting at Q. 
Fluctuations induced by these nearby phases could be captured within an RPA treatment 
of a model with a more realistic bandstructure. 

For simplicity we have discussed the two channel Anderson model involving fluctuations 
from a 5/ 2 T 5 ground state to 5/ 1 (J = 5/2). However, the more realistic case involves 
fluctuations to 5/ 3 , whose low energy states have J = 9/2, and are split into five Kramers 
doublets by the tetragonal crystal field, 

\T\±) = a\ ±5/2) +b\ =f3/2) 
\T 2 7 ±) = -b\ ±5/2) +a\ =F3/2) 
|rj' 2 ' 3 ±> = c 1 ' 2 ' 3 ! ± 9/2) + d 1 ' 2 ' 3 ! ± 1/2) + e 1 ' 2 ' 3 ! =F 7/2). (37) 

There are two generic situations: either a T 7 doublet will be lowest in energy, and the valence 
fluctuations will be determined by the overlap, 

|r 7 ±) = a4 T \r 5 ±) + M T |r 5 ±), (38) 

or a T 6 doublet will be lowest in energy, with the relevant overlap, 

|r 6 ±) = a4^\r 5 ±) + (34 t \T 5 t), (39) 

where the form-factors are as above. In both cases fluctuations will involve conduction 
electrons in both T 6 and T 7 symmetries. When the lowest excited state is a r 7 , the valence 
fluctuation Hamiltonian is given by, 

HvfzU) = v^j r6T |r 5 ±)(r 7 ± | + v^j r7T |r 5T )(r 7 ± | + h.c. (40) 
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Now we can follow the same slave boson procedure as discussed previous, where now 
represents a 5/ 3 state and carries a negative charge, and the valence fluctuation term is the 
particle-hole conjugate of the 5/ 1 case, cj ccr V crQ (k, j)XaO) e_ * k Rj - ^ we perform a particle- 
hole transformation on x, we regain the Hamiltonian, (34) with only a sign change for A k . 



V. MAGNETIZATION AND SUSCEPTIBILITY CALCULATION 

The uniform basal plane conduction electron magnetic susceptibility is 

X xy = ~(9^b) 2 T J2J2 Ti \° X V C ^ k + Q,ioj n )a y G c (k + Q, k, iu n )\ , (41) 
while the staggered conduction electron magnetization is 

m c (Q) = - gi i B T E Tr ^ C ( k ' k + Q- ™n) (42) 

iu> n k 

In order to calculate these, we require the full conduction electron Green's function, which 
can be found from the Hamiltonian by integrating out the f-electrons, 

, / ioo n — e k \ / iu n — A k \ . 

y iu n - e k+Q J \ lUJn ~ Ak +Q / 

where 

V k =f V « k M. (44) 
\-Vm -Vat J 

Using isospin, f to represent k, k+Q space, we split the conduction electron energy, e k into 
e ok — \ ( e k + e k+Q)j tik — \ (e k — e k+ Q) into the particle-hole symmetric and antisymmetric 
parts (and similarly with A ok , A ik ). So now we can write the conduction electron Green's 
function as: 



[£ c (k, = (iu n - e ok ) - e lk r 3 - V p^ A x ° k ^ Alk x ? Vt (45) 

(lL} n — A ok J — A lk 



For completeness, the f-electron Green's function will be, 



[0'(k, iu)] 1 = (iu n - A ok ) - A lk r 3 - V ^ n ,2 lk ? V. (46) 
L J - e 0k ) 2 - ef k 
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The Green's function can be expanded in terms of r <g> a, where a describes the spin. To do 
so, we take TrVVV a er 6 and TrVT 3 VV a <7&. The relevant parameters are: 



V£ + = TrV k Vj = 2Tr 
V k 2 _ = TrV k r 3 Vjr 3 = 2Tr 
A k _ = Tr V k r 3 V k r 2 = 2iTr 
A k+ = TrVkVtna = 2Tr 



VekVi + V 7k V} k 



VeicV^ - V 7k V. 



t 



7k V 7k 



v 6k v 7k - v 7k vl 
v 6k v} k + V 7k V 6 f k 



a 



(47) 



All other combinations are uniformly zero. The conduction electron Green's function can be 
written in the form [G ] -1 = Ar + £?t 3 + C ■ Bt\ + Dr 2 . The eigenvalues are found by taking 
det [{? c (k, u))] 1 = 0, leading to an eighth order polynomial. For the special particle-hole 
symmetric case where eo = Ao = 0, this reduces to u A — 2a k w 2 — 7 k , which can be solved 
analytically. The four (doubly degenerate) eigenvalues are E kn and are found numerically 
on a grid of k points. Due to the structure of the Green's function, we can write, 

1 



G c (iu n , k 



Atq — Bt 3 — C ■ (JT i — Dr 2 



Yl v (iu) n - E kr) ) 
A = (iu n - e ok ) [{iuJn ~ A 0k ) 2 - A 2 k ] - (iu n - A 0k )F k+ 

B = -ei k [(iu n - A 0k ) 2 - A 2 k ] - Ai k \/ k 2 _ 

C = —{iu n - A ok )A k+ 

D = -A lk A k _ 

which makes it particularly easy to write down the conduction electron magnetization, 



(48) 



m, 



: (Q) = -(^ B )T^Tr[^ c (k,^ n )r 1 ] 

— — (p/x ) ^ kr; ~ ^ ok ^ 



f(E kT) )A k+ 



(49) 



and susceptibility, 

X xy = -(g^ B ) 2 Tj2J2^^ gC ^ iu >y gC ^ iuj -)} 

iw n k 

2(-E kr? - A k)/(-E'kj?) + (-Ekj? - Ao k ) 2 /'(-^k»)) 



-(<W*) 2 £ 

kr; 

-E 



2(-E'k»7 — Aok) 2 /(-Ek»7) 



A x A y 



(50) 



Note that above integral may be positive or negative. The functions / and /' are the Fermi 
function, f(x) = (e~ x l T + l) 1 and it's derivative, f'(x) = df(x)/dx, respectively. 
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Although the conduction electrons develop a magnetic moment, the f-electrons have no 
dipolar or quadrupolar moments, rhf = 0. Upon closer examination, the quadrupolar mo- 
ments are found to vanish because of a d-wave form factor, implying that there is no net 
quadrupole moment and thus no associated lattice distortion, even for uniform hastatic or- 
der. As a d-wave quadrupole is a hexadecapole, this ultimately means that like Haule and 
Kotliar[5], we have staggered (J x J y + J y J x )(J^ — Jy) hexadecapolar moments. However, 
unlike Haule and Kotliar, where the hexadecapolar moments are the primary order param- 
eter (and thus of order one), here the hexadecapolar moments are a seconday effect of the 
composite hastatic order, and like the conduction electron moments, will be of order T K /D. 
The f-electron moment has an identical form to the conduction electron moment, (37), with 
e <H> A everywhere, and the relevant form-factor becomes A£ + = TrVkoViV^. Given how 
difficult it is to observe large hexadecapolar moments, the hexadecapolar moments associ- 
ated with hastatic order will almost certainly be unobservably small. By contrast, in the 
antiferromagnetic phase, the f-electrons develop a large c-axis magnetic moment. 



VI. CURRENT NEUTRON SCATTERING BOUNDS ON THE TRANSVERSE 
MOMENT 

We predict a staggered basal plane magnetic moment on the order of .01/ie/U. However, 
such a moment has not been previously detected, despite a large number of neutron studies. 
Early neutron studies focused on the small, .02 — .04/ie c-axis moment [13, 14], which was 
later shown to be due to small, local regions of large moment antiferromagnetism[15]. Such 
inhomogeneous antiferromagnetism persists even in high quality single crystals, making it 
essential to identify the moment orientation. The moment orientation has not been examined 
since these early works[13, 14], whose resolution of .01/ie is on the order of our prediction. 
As the neutron signal is proportional to the component of the moment perpendicular to Q: 
m_|_ oc Q x (m x Q), measuring along Q = (001) will isolate the transverse magnetic moment. 
More recent studies have focused on Q = (100), where the extrinsic c-axis moment is always 
seen[12]. Detection of our predicted mi will require high resolution, likely spin-polarized, 
measurements along Q = (001). 
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VII. DENSITY OF STATES 



To simplify the calculation of the ground-state properties, we have taken e k = — A* an d 
Aok = A to be constant, and chosen the filling of the conduction sea so that at T = 0, 
A + fj, = 0. In this case, the Hamiltonian %(k) Qj g in (34) can be diagonalized analytically, 
leading to four doublets \krjcr) (a = ±1), with eigenvectors 

u a (kr]cr) = (a\kr)a) (51) 

s where % Qj g(k) •-u(k^o-) = E kri u(kr]a),. The degeneracy of the eigenvalues is guaranteed by 
the invariance of the Hamiltonian under time-reversal plus a translation by Q. The energies 
E kv are given by 



£^ = ^±^±^-^1 (52) 

where a k = V£ + + \ (e? k + A? k ) and 7 ^ = (e lk A lk - V*_) 2 + (A k _) 2 . 
The total density of states is then given by 

A M = ft,) = -Im y £ (53, 

k?j v y »7 ' 

where the integral is over the Brillouin zone. Numerically, this quantity was computed by a 
sum over a discrete set of momenta, dividing the Brillouin zone into 40 3 points and using a 
small value of 5 to broaden the delta-function into a Lorentzian. 



VIII. G-FACTORS 



The Zeeman energy is determined by the Hamiltonian 

-B-M = - ^lM^-B 



(54) 



where Vk = (c k , c k+Q , Xk, Xk+q) T and 



M = 



(2fi B a \ 

1\l b o 

g fl i B a z 

V g f HBO*) 



(55) 
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where gf is the effective g-factor of the Ising Kramers doublet. In a field, the doubly- 
degenerate energies, |k?7cr) (a = ±1) are split apart so that AE kv = \E kv ^ — E kv i\ = g kv (9)B, 
so the g-factor is given by g^O) = dA J^ v ■ Now we are interested in the Fermi surface 

B — >0 

average of the g-factor, given by 

9{e) = ^jm 

These quantities were calculated numerically, on a 40 3 grid, using gf = 2.9 for the effective g- 
factor of the local Kramers doublet. The resulting g-factor in the z— direction is reduced to 
9eff(@ = 0) = 2.6 because of the admixture with conduction electrons. The delta-functions 
were treated as narrow Lorentzians 6(E) = \lm(E — iff)" 1 , where r\ is a small postive 
number. The g-factors at each point in momentum space were computed by introducing a 
small field 6B into the Hamiltonian, with the approximation g^O) = \E^ — E^\/SB. 



IX. MODEL TUNNELING CONDUCTANCE AND NEMATICITY CALCULA- 
TION 

To calculate the tunneling density of states, we assume that the differential conductance 
is proportional to the local Green's function on the surface of the material 

^(x)ocA(x,eV) (56) 

where 

-| poo 



1 f°° 

A(x,u) = -ImG aa (x,u-i6)=J2 / ^({^(x,t),^(x,0)}> 

77 „ J — oo 



e 



is the imaginary part of the local electronic Green's function. To calculate this quantity, we 
decompose the local electron field in terms of the low energy fermion modes of the system, 
writing 

Vv(x) = M x " R i\) c i° + <^*( x - R i)fjre« + " R-;)/ir 7 «] (57) 

j 

where c (|x — R,-|) is the wavefunction of the conduction electron centered at site j, while 

0L = 6 (|x-R,|)^L(x-R,) 

0L = 7 (|x-R,|)^ 7 a (x-R,) (58) 
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are the wave functions of the T 6 and T 7 - f-orbitals centered at site j. 

Projected into the low energy subspace, we have fjr 6 a — > ({Bj)Xj)a an d fjr 7 a — > 
((Bj)aiXj)a- Writing £>j = bUj, and Xj — UjXj these expressions become fjv 6 a - > bxj 
and fjr 7 a b(n ■ (T)e _lQ ' Rj Xj. Making these substitutions and rewriting the field operators 
in momentum space we obtain 



\ 



C k+Qct 
VXk+Qa/ 



(59) 



where 



A,(x) = ( c (|x|)<W _ * Q ' R ^c(|x|)<W 

b€a(x) + e-^' R ^0L(x)(n^),e-^- R ^0L(x) + (60) 

We choose a layer where e -1 ^' 1 *' 3 **' = +1, then on this layer the local Green's function is 
given by 

(61) 



where 



G(x, u) = ^(x - Rj) ■ Qji{u) • At(x - R,) 

3,1 

Q 5 M= Tr[(l + r 1 )^(k,o;)] e - jk -( R ^ R ) 

is a trace only over the momentum degrees of freedom, so Qji is a four by four matrix for 
each pair of lattice points j and /, where A(x) = (<f) c (\x.\)5 aa , 60^ a (x) + 60£ a (x)(n ■ a)). 
The final spectral function is then 



A(x.,u) = -Im Tr 



7T 



A(x - Rj) ■ Gji(u - iS) • A f (x - Rj) 



(62) 



To evaluate this quantity, the summations were limited to the four nearest neighbor sites at 
the corner of a plaquette. The positions x were taken to lie in the plane of the U atoms. 
The wavefunctions 6 (|x|) = e~l x l/ a , 7 (|x|) = e l^l/* 3 - and c (|x|) = e~^^ a were each taken 
to be simple exponentials of characteristic range equal to the U — U spacing a. 

The nematicity of the tunneling conductance was then calculated numerically from the 
spatial integral 

f A(x, eV) sga(xy)dxdy 



V (eV) = 



(7 dxdyA{x, eV) 2 - [/ dxdyA{x, eV)] 2 ) 



1/2 



(63) 
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